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SUMMARY 

The  response  of  any  Optical  System,  used  in  either  a  weapon  or  surveillance  mode,  is 
dependent  upon  the  emission  of  electromagnetic  energy  from  born  target  and  background 
occupying  the  field  of  view.  The  background  not  only  contributes  to  the  infrared  scene 
(clutter)  but,  additionally,  by  reflection  and  absorption,  can  significantly  contribute  to 
the  target's  signature. 

Since  the  maritime  environment  can  be  conceived  as  containing  an  opaque  homogeneous 
body  of  water,  its  signature  will  depend  upon  the  predictable  behaviour  of  such  a 
surface  in  response  to  the  natural  elements. 

This  report  proposes  a  statistical-numerical  method  which  can  be  used  to  obtain  the 
infrared  signature  of  both  a  calm  and  wind  roughened  sea  in  the  1  to  15  gm  waveband. 
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1.  INTRODUCTION 

To  determine  the  response  of  any  optical  system  (weapon  or  surveillance)  to  the  electromagnetic  emission 
of  a  target  under  consideration,  it  is  often  necessary  to  have  knowledge  of  the  target's  background.  The 
background  can  not  only  contribute  directly  to  the  scene  (clutter)  but  additionally  by  interaction  with  the 
target  it  can  contribute  signficantly  to  the  target's  signature. 

The  maritime  environment,  as  a  background  scene,  can  be  conceived  as  a  homogeneous  body  of  water  that 
responds  to  the  natural  elements  in  a  predictable  manner.  Because  of  such  behaviour  the  maritime 
environment  conveniently  offers  itself  to  be  mathematically  modelled  such  that  within  the  1  to  14  pm 
waveband  its  infrared  signature  can  be  determined. 

Two  approaches  can  be  used  to  model  the  maritime  environment.  The  first  is  a  statistical  model  which 
was  developed  by  Saunders(ref.l)  to  predict  the  radiance  of  the  sea  in  the  12  to  14  pm  waveband.  The 
second  approach  is  a  method  used  to  predict  the  sea's  infrared  signature  by  a  geometrical  interpretation  of 
the  sea  surface.  This  method  offers  both  a  spatial  and  time  variant  infrared  scene  of  a  calm  and  wind 
roughened  sea.  Although  this  approach  may  seem  to  have  a  definite  advantage  over  the;  of  a  static 
statistical  surface,  if  the  radiance  of  the  sea  is  needed  over  a  large  portion  of  the  surface  (many  wa  v  e 
periods)  then  the  computing  time  necessary  to  evaluate  such,  can  become  both  timely  and  costly.  This 
situation  can  arise  when  an  imager  for  example  contains  a  wide  field  of  view  or  when  a  narrow  field  of 
view  is  used  at  an  extended  range  from  a  region  of  the  sea  surface.  Such  a  model  will  be  discussed  in 
another  report(ref.2). 

The  following  report  discusses  the  statistical  approach  used  to  predict  the  signature  of  both  a  calm  and 
wind  roughened  sea.  It  combines  the  method  of  Cox  and  Munk(ref.3)  and  that  of  Saunders  to  form  the  basis 
onto  which  a  numerical  method  is  developed,  to  produce  a  signature  over  the  entire  1  to  15  pm  waveband. 


2.  RADIANCE  MODEL 

The  signature  of  the  sea  surface  in  the  8  to  12  pm  waveband  is  well  documented  and  has  been  treated  by 
several  authors.  The  requirement  is  to  supplement  and  expand  on  such  models  to  include  the  shorter 
wavebands  so  that  the  usable  range  extends  from  1  to  14  pm.  The  difficulty  encountered  using  shorter 
wavelength  regions,  is  that  the  reflected  radiation  from  the  sea  surface  which  originates  from  the  sky 
significantly  becomes  a  function  of  azimuth  in  addition  to  elevation,  whereas  the  sky  radiance  in  the 
longer  8  to  12  pm  waveband  is  usually  treated  as  only  a  function  of  elevation.  The  azimuth  dependence 
must  be  accounted  for  and  is  detailed  in  the  following  model  description. 

2.1  Calm  sea 

The  statistical  model  closely  follows  the  approach  taken  by  Saunders  such  that  if  (B)  represents  the 
Black  Body  radiance  of  the  sea  and  (NS)  represents  the  radiance  of  the  sky,  then  at  observation 
angle  <b  (see  figure  1)  the  total  radiance  of  a  calm  sea  is  given  by 


N  =  (B-R(B-NS)) 


(1) 


a  =  Azimuth  of  Max  Slops 
3  =  Angls  of  Max  Slops 

<p  =  Obssrvsr  Angls  to  Y-Axis 
H  =  Reflected  Sight  Lins  Zsnith  Angls 
v  =  Rsflsctsd  Sight  Lina  Azimuth  Angls 
2<o  =  Angle  between  Observer  and  Reflected  Sight  Line 

Figure  1.  Geometry  for  reflection  angles  from  a  sea  surface  facet 
where  R  is  the  reflectance  for  the  sea  and  is  calculated  from  the  standard  Fresnel  equation 

2  R(to)  =  sin2  ((o'  -  e)  esc2  (o’  +  e')  +  tan2  (co'  -  e')  cot2  (o'  +  e’) 

where  cu'  and  e'  are  the  incidence  and  transmittance  angles  respectively. 

2.2  Wind  roughened  sea 

In  the  presence  of  wind,  equation  (1)  has  to  be  modified  to  account  for  the  waves  of  the  sea  which,  for  a 
given  look  down  angle'  (<)>),  will  reflect  portions  of  the  sky  extending  over  the  entire  hemisphere. 
That  is,  if  Q  =  0  corresponds  to  looking  at  the  horizon,  then  for  0  =  80°  (that  portion  of  sea  below  the 
observer)  there  is  a  finite  probability  that  in  the  presence  of  waves  there  will  be  a  contriubtion  of  sky 
radiance  from  the  horizon  albeit  small.  Because  a  probability  or  probability  distribution  is  used,  the 
resulting  sea  radiance  is  a  time  average  of  the  instantaneous  radiance. 

The  wave  slope-distribution  function  p(Zx,  Zy)(ref.3)  is  defined  as  the  probability  of  slope 
components  Zx,  Zy  occurring  in  a  unit  area  of  horizontal  sea  surface.  Thus  p(Zx,  Zy)  6Zx  5Zy  represents 
the  fraction  of  surface  with  slope  components  in  the  interval  Zx±8Zx/2  and  ZythZy/i  and 
equation  (1)  becomes 
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N  =  cosec  <pf  f  (B  -  R  [B  -  NSJ)  cos  w  see  (3  p(Zx,  Zy)  dZx  dZy 


(2) 


where  by  referring  to  figure  1  and  Appendix  I,  and  for  a  look  down  angle  0,  a  (3  co  p  and  v  can  be 
determined.  Zx  and  Zy  are  the  slope  components  in  the  X  and  Y  directions  respectively.  The  term 
coscec  <(>  is  included  as  a  constant  for  a  given  look  down  angle  on  the  right  hand  side  of  equation  (2)  so 
that  N  becomes  the  radiance  with  respect  to  line  of  sight. 

To  account  for  the  dependence  of  radiance  on  the  azimuth  direction  (v),  equation  (2)  becomes 


N(<J>)  =  Scosecp 


tan 


I  J 


{B  -  R  [B  -  NS(p,v)])  cos  co  Sec  (3  p(Zx,  Zy)  dZx  dZy 


(3) 


The  factor  (S)  in  equation  (3)  is  used  to  account  for  the  reduction  in  the  radiance  due  to  the  effects  of 
both  shadowing  and  multiple  reflections.  The  finite  integration  limit,  tan  0,  accounts  for  including 
only  those  slope  components  which  do  not  exceed  the  inclination  of  the  line  of  sight. 

2.2.1  Shadowing 

The  derivation  and  discussions  of  the  slope-shadowing  factor  can  be  found  in  detail 
elsewhere(ref.l),  accordingly  it  is  sufficient  to  only  provide  the  expression  for  S  which  accounts 
for  die  fraction  of  the  surface  which  remains  hidden  from  the  line  of  sight  due  to  the  masking  of 
waves  by  other  waves  closer  to  the  observer.  S  is  related  to  the  error  function  and  is  given  by 


S  =  2  (1  +  erf(V)  +  (VVrc)'1  exp  -  V2}'1 


where  V  =  8"1  tan  $  and  8  is  the  root  mean  square  slope,  the  value  which  is  obtained  by  the 
expression(ref.3) 


S2  =  0.003  + 5.12  x1CT3W 


W  is  the  wind  velocity  in  ms"’.  To  a  first  order  the  slopes  of  the  sea  surface  are  normally 
distributed  and  independent  of  wind  direction  thus  the  wave  slope  distribution  function  is  related 
to  the  wind  speed  by  the  following  expression 
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p(Zx,  Zy) 


The  statistical  sea  radiance  equation  (3)  does  depend  on  wind  speed  but  not  wind  direction. 

2.22  Multiple  scattering 

Multiple  scattering  arises  when  the  ray  of  electromagnetic  radiation  travelling  towards  the 
observer  has  undergone  at  least  two  reflections.  This  condition  occurs  when  the  elevation  angle  p 
in  figure  1  has  the  property 


p  >  90° 


that  is  the  ray  reflected  toward  the  observer  does  so  by  a  reflection  such  that  the  incident  ray 
originates  from  below  the  horizon. 

Multiple  scattering  is  handled  iteratively  in  the  following  report  by  a  "Successive 
approximation”  method  which  can  be  described  by  the  following.  Sea  radiance  values  are  first 
calculated  for  both  normal  and  grazing  viewing  angles.  This  provides  two  boundary  values  such 
that  the  contribution  of  radiance  from  multiple  scattering  is  calculated  initially  by  interpolation 
between  these  boundaries  resulting  in  the  determination  of  the  total  sea  radiance  from  a  point  just 
below  the  grazing  viewing  angle  of  the  horizon.  Having  now  three  points  to  the  sea  radiance,  the 
next  radiance  value  can  be  determined  by  appropriately  linearly  interpolating.  In  such  a  manner 
each  point  of  radiance  is  determined  for  angles  progressively  lower  than  the  grazing  viewing 
angle  and  toward  the  normal  viewing  angle.  Thus  the  iterative  technique  builds  up  the  sea 
radiance  map,  starting  from  the  horizon  and  progressively  determining  points  of  radiance  towards 
the  normal  viewing  angle. 

The  effect  of  using  a  boundary  at  grazing  viewing  angle  which  was  either  too  large  or  too  small  in 
calculating  the  second  sea  radiance  value,  was  checked  by  deliberately  using  the  largest  possible 
radiance  value  and  recording  the  sea  radiance,  and  then  the  smallest  possible  radiance  vlaue  (the 
Zenith  radiance)  anri  recording  the  value.  The  difference  was  found  to  be  less  than  1%. 


3.  INTEGRAL  ELEVATION 

Equation  (3)  cannot  be  easily  analytically  integrated.  It  was  therefore  necessary  to  numerically  evaluate 
the  integrals.  With  the  appropriate  change  of  variables  (Appendix  ID  equation  (3)  reduces  to  the  form 


N  (0,0)  =  k(I,  +  y 


(4) 
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where 


oo 

h  =  82  j 

-oo 


oo 

J  G  (y,y)  e-V  dy 
o 


e"?  2  dy 


(5) 


The  inner  integral  can  be  evaluated  using  Gauss-Laguerre  Quadrature  and  if  we  denote  the  evaluation  of 
the  inner  integral  as  H'(y)  then  equation  (5)  reduces  to 


OO 

I,  =  52  J  H’(y)e-Y 2  dy 

-oo 


(6) 


which  is  in  a  form  suitable  to  be  evaluated  by  Gauss-Hermite  Quadrature. 

By  a  further  manipulation  in  terms,  the  second  integral  in  equation  (4)  can  be  reduced  to 


h 


oo 


52  J  e 


1 

/  F(5y,6(e  + pe))e-£2('+P)2dp 

-1 


e"T 2  dy 


Again  the  inner  integral  can  be  evaluated  using  Gauss-Legendre  Quadrature  and  representing  the 
evaluation  by  M(y),  I2  reduces  to 


00 

I2  =  52  j  M(y)  e"Y 2  dy 

-00 


which  is  of  the  form  of  equation  (6)  and  can  be  similarly  evaluated. 

As  discussed  earlier,  to  evaluate  equation  (4)  to  account  for  any  waveband  then  rather  than  evaluating 
the  sky  radiance  NS  (p,v)  semi  empirically,  the  radiance  information  is  read  from  a  data  file  which 
contains  values  for  the  whole  hemispherical  sky.  Values  which  are  needed  between  the  data  file  data 
points,  are  found  by  interpolation  of  the  data  set.  At  present  it  has  been  sufficient  to  use  only  linear 
interpolation. 

The  additional  advantage  of  using  such  a  file-reading  technique  is  that  experimentally  measured  data 
sets  of  sky  radiance  values  can  be  entered  if  they  are  available.  A  flow  diagram  is  shown  in  figure  2  to 
obtain  the  total  sea  radiance  given  wind  velocity,  azimuth  of  the  observer  and  the  sky-radiance  data 
file. 


i 


Figure  2(a).  Row  diagram  of  radiance  model  mam  program 
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SUBROUTINE  ROUGH 


Figure  2(b).  Row  diagram  of  subroutine  rough 


ERL-0473-TR 


8 


4.  RESULTS  AND  DISCUSSIONS 


4.1  Long  waveband 

To  confirm  the  validity  of  the  numerical  technique  used  to  read,  interpolate  and  integrate  within  the 
model,  comparison  of  the  sea  radiance  from  the  model  oceanrad  was  made  in  the  8.3  to  12.5  pm 
waveband  against  the  values  obtained  by  Saunders  as  shown  in  figures  3(a)  and  3(b).  Similarly  in 
this  waveband,  some  important  conclusions  must  be  comn  tented  on: 

(i)  As  the  'look-down'  angle  (shown  as  negative  angles)  decreases  from  90°  to  45°,  the  sea 
becomes  only  fractionally  cooler  with  only  a  small  departure  from  the  Black  Body  temperature  of 
the  sea.  Even  at  90°  the  sea  radiance  will  always  be  smaller  than  that  of  a  Black  body  since  its 
reflectivity  is  always  greater  than  zero. 

(ii)  Except  for  a  condition  of  zero  wind  speed,  the  sea  radiance  decreases  with  the  magnitude  of 
the  look  down  angle.  The  radiance  at  the  horizon  for  each  value  of  wind  speed  varies 
discontinuously  to  that  of  the  sky  above.  The  magnitude  of  discontinuity  decreases  for  increasing 
wind  speed  hence  sea  roughness. 

(iii)  The  calm  sea  exhibits  no  discontinuities  at  or  near  the  horizon,  however,  the  sea  radiance, 
within  an  angle  of  0.5°  of  the  horizon,  becomes  increasingly  hotter  until  it  reaches  the  near 
ambient  temperature  of  the  horizon  sky.  Because  the  calm  sea  is  flat,  then  for  grazing  look  down 
angles,  the  warm  horizon  sky  will  be  reflected  to  the  observer.  However,  in  the  presence  of  waves 
and  as  a  result  of  the  Shadowing  Factor  S  in  equation  (3),  then  the  sea  appears  dish  shaped'. 
That  is  for  grazing  look  down  angles  the  sea  now  has  a  non-zero  slope  probability  p(Zx,  Zy)  and 
the  colder  sky  is  reflected  accordingly. 

Figure  4  shows  the  mean  radiance  of  the  sea  in  the  8  to  12  pm  waveband  as  a  function  of  the  refractive 
index.  Because  different  authors  have  used  different  refractive  indices  for  sea  water  (Cox  and 
Munk  1.338,  Saunders  1.237  and  others  which  obtain  a  range  of  values  between  1.334  to  1.36  depending 
upon  the  salinity  and  chlorinity  content)  a  large  variation  of  the  sea  radiance  may  result  from 
differenting  indices.  The  three  sets  of  curves  in  figures  4(a)  and  4(b)  represent  the  effect  of  wind 
speeds  1.36  ms"1  and  17  ms"1  which  correspond  to  sea  roughness  of  8  =  0.1  and  0.3  respectively.  Most 
variation  in  radiance  occurs  at  grazing  look  down  angles  for  both  maximum  and  minimum  velocities. 
Table  1  contains  the  numerical  value's  of  figures  4(a)  and  4(b)  for  the  largest  and  smallest  values  of 
refractive  indices.  The  diff<'  ence  in  r'diance  for  such  values  can  be  seen  at  the  bottom  of  the  table 
indicating  that  an  arbit."  hoice  of  refractive  index  within  the  given  range  will  result  in  a 
maximum  variation  in  r  ;  ee  of  less  than  5%.  This  percentage  will  decrease  as  the  magnitude  of 
the  look  down  angle  increases  is  therefore  expected  that  an  incorrect  choice  of  the  refractive  index 
within  the  specified  ra  ige  wii<  ,.ot  significantly  change  the  sea  radiance. 

4.2  Short  waveband 

As  discussed  earlier,  the  sky  radiance  in  the  short  infrared  waveband  is  a  function  of  azimuth  angle. 
The  main  contribution  of  angular  dependence  comes  from  direct  and  indirect  (scattered)  solar 
radiation.  As  shown  in  figure5,  elevation  is  plotted  against  mean  sky  radiance  for  five  different 
azimuti.  angles. 
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REFRACTIVE  INDEX  (AIR/SEA) 

Figure  4(a).  Mean  radiance  against  refractive  index  for  8  =  0.1 


REFRACTIVE  INDEX  (AIR/SEA) 


Figure  4(b).  Mean  radiance  against  refractive  index  for  8  =  0.3 


VIEWING  ANGLE  (DEGREES) 
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Such  sky  radiance  onJy  approximates  the  expected  3  to  4  pm  radiance  map,  since  it  is  the  purpose  of 
this  report  to  present  a  model  and  not  to  present  a  comparative  assessment  between  measured  and 
numerically  calculated  infrared  radiance  maps.  As  shown  by  the  increase  in  radiance,  the  sun  is 
approximately  43°  in  the  east. 

Figure  6(a)  shows  both  the  sea  and  sky  radiance  as  a  function  of  'look  down  angle’  for  an  azimuth 
direction  of  180°  (directly  away  from  the  sun)  and  figure  6(b)  shows  the  radiance  for  an  azimuth 
direction  of  0°. 

In  each  figure  the  two  curves  represent  the  calculated  radiance  for  a  wind  velocity  of  0  and  20  ms'5. 
The  radiance  curve  in  figure  6(b)  for  0  ms"*  is  similar  to  the  LWB  results  'in  trend'  in  that  the  values 
decrease  as  depression  decreases  to  a  value  at  -5°  then  rises  sharply  at  the  horizon.  However,  it  does 
contain  the  additional  radiance  of  scattered  sun  radiance  which  can  be  seen  reflected  in  the  sea  at  -43°. 
The  intensity  from  such  a  reflection  is  reduced  as  expected  from  the  increased  emissivity  of  the  sea  for 
that  angle. 

For  sea  state  or  roughness  other  than  6  =  0,  the  short  waveband  radiance-profile  can  differ 
significantly  from  that  of  the  long  waveband.  As  shown  in  figure  6(b)  (depending  upon  sun  elevation), 
a  situation  can  be  attained  where  the  radiance  remains  constant  or  may  even  increase  toward  the 
horizon.  This  results  in  the  discontinuity  fo  the  sea/sky  interface  either  not  being  present  or  inversed 
to  that  which  is  usually  expected. 

Figure  7  represents  radiance  profiles  for  wind  velocities  of  20  ms'1  for  a  number  of  different  azimuth 
angles.  This  demonstrates  the  totally  different  sea  signatures  which  can  be  achieved  depending  upon 
the  observed  azimuth  angles.  Notice  that  for  look  down  angles  with  a  magnitude  greater  than  45°  the 
radiance  of  the  sea  approaches  that  of  a  'black  body’. 

43  Image  representation 

The  data  output  from  the  sea  radiance  model  was  written  into  a  data  file  of  37  columns  and  64  rows. 
Using  two  dimensional  linear  interpolation,  a  second  data  file  of  128  columns  and  64  rows  was 
generated.  In  this  format  the  radiance  information  could  be  coded  and  displayed  on  a  PC  (16  bit 
Amiga)  using  an  identical  'grey  scale'  display  as  that  which  the  Aga  Thermavision  Termal  Imager 
formats.  At  present  the  coding  is  such  that  the  grey  scale  is  limited  to  11  levels  for  display  purposes. 

For  a  single  azimuth  angle,  the  radiance  program  takes  10  s  to  run  for  a  zero  input  wind  speed,  and  90  s 
for  any  finite  wind  speed  other  than  zero.  Numerically  generated  radiance  images  of  the  maritime 
scene  using  the  sea  radiance  program  and  the  long  and  short  waveband  (LWB  and  SWB)  are  shown  in 
figures  8  and  9.  Each  image  results  from  the  observer  subtending  a  60°  field  of  view.  The  horizon  in 
each  case  occurs  at  the  midpoint  of  the  vertical  field  of  view. 

In  particular,  figure  8(a)  shows  the  8  to  12  pm  (LWB)  image  for  a  wind  speed  of  zero  (left  hand 
image),  and  5  ms"1  (right  hand  image).  In  each  picture,  and  to  the  extreme  right,  a  grey  scale  of 
11  levels  is  displayed.  'Hotter  regions'  are  characterised  by  a  'lighter'  level  of  grey  scale.  The 
horizon  separating  the  colder  sky  from  the  warmer  sea  can  be  clearly  distinguished. 
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Figure  7.  SWB  sea/sky  radiance  profile  for  different  azimuth  angles 

Figure  8(b)  shows  the  numerical  model  of  a  Guided  Missile  Destroyer  inserted  onto  the  LWB  sea 
background.  The  grey  scale  of  the  target  at  this  stage  is  not  meant  to  be  compatable  to  its  background 
since  such  an  image  is  included  to  only  demonstrate  an  example  of  an  application  of  the  sea  radiance 
imaging  facility. 

The  left  picture  of  figure  9(a)  represents  the  8  to  12  pjn  (LWB)  image  which  can  be  compared  to  the 
right  picture  which  represents  the  SWB  image  of  the  maritime  scene  for  a  common  wind  velocity  of 
zero.  The  azimuth  dependence  of  the  radiance  can  be  clearly  seen  for  the  SWB.  Figure  9(b)  shows  the 
effects  of  different  wind  speeds  of  zero  and  5  ms*1  for  the  SWB  image.  Again  variation  in  radiance  is 
evident  as  a  function  of  both  wind  speed  and  azimuth  direction. 

4.4  Irradiance 

Once  the  radiance  of  the  sea  surface  has  been  determined,  the  incident  power  of  the  electromagnetic 
energy  (irradiance)  at  a  unit  surface  on  the  observer's  optics  can  be  determined.  If  the  range  between 
the  emitting  sea  and  observer  is  sufficiently  large  (at  wavelengths  corresponding  to  the  atmospheric 
windows)  or  if  wavelengths  within  the  absorption  profiles  of  the  atmosphere  are  considered  then 
effects  from  the  intervening  atmosphere  must  be  account  for.  This  is  considered  in  the  following 
irradiance  equation  which  is  given  by: 


HO$)  =  ^  [N(4>>  t(X>  +  Na  (0,R)1 


(7) 
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where  N(0)  corresponds  to  the  sea  radiance  as  given  by  equation  (3)  and  Na(0,R)  is  the  total  radiance 
of  the  intervening  atmosphere.  It  must  also  be  noted  that  the  irradiance  is  more  accurately  the 
spectral  irradiance,  although  in  the  atmospheric  windows,  t(X)  the  atmospheric  transmittance,  is 
very  nearly  constant.  However  if  t  is  a  function  of  the  wavelength  and  appreciably  changes  over  the 
waveband  of  detectability,  then  the  observed  irradiance  is  given  by 


H"»  -  s 


\2 

j  D(X)  dX 


where  AX  =  X2  -  Xi  and  represents  the  waveband  in  which  the  observation  is  made.  The  term  D(X) 
accounts  for  the  transmittance  for  the  spectral  passband  of  the  detector. 

The  term  is  parenthesis  in  equation  (7)  differs  from  the  sea  radiance  N(<)>)  with  a  magnitude  which 
depends  upon  both  the  atmospheric  transmittance  t(X)  and  the  airlight  emission  radiance  Na(X). 
Therefore,  in  the  3  to  4pm  and  8  to  12pm  waveband,  effective  sea  temperatures  can  differ 
significantly  from  the  actual  combined  sea/sky  temperatures  for  the  same  look  down  angle.  This  may 
well  be  of  significance  for  a  weapon  imaging  system  tracking  a  target  at  or  near  die  horizon. 

Figure  10  shows  an  optical  surviellance  system  relative  to  the  horizon,  where  r  is  the  earth  s  radius, 
s  is  the  vertical  hieght  of  the  observer,  Lm  is  the  distance  from  the  observer  to  the  horizon  and  L2  is 
the  line  of  sight  to  the  sea  surface  which  makes  an  angle  p  to  the  vertical.  From  the  simple  geometry, 
the  angle  from  the  vertical  to  Lm  is  given  by: 


0  =  sin-1 


hence 


Lm  =  r  cot  $ 


and  using  the  cosine  law  for  triangles,  and  simplifying,  L2  is  given  by: 


L2  =  (r  +  s)cosp- {r2-sin2p(r  +  s)2}1/2 
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Figure  11  shows  a  plot  of  L2  as  a  function  of  viewing  angle  p,  for  altitudes  (s)  of  10  m  and  1000  m 
respectively.  For  an  observers  at  10  m,  a  look  down  angle  to  the  horizon  gives  an  L2  of  272  m.  The 
average  atmospheric  transmission  for  this  range  is  0.91 1  so  that  the  horizon  contrast  will  still  be 
clearly  distinguished.  However,  at  an  altitude  of  1  km,  even  though  the  decrease  in  path  length  as  a 
function  of  look  down  angle,  is  substantial,  for  p  =  75.5°  L2  =  4.0  km  which  corresponds  to  a 
transmissivity  of  only  0.25.  Therefore,  depending  upon  the  atmospheric  conditions  and  observers 
height  above  sea  level  there  may  be  very  little  or  no  discemable  contrast  between  the  sky  and  sea, 
furthermore  the  sea  may  appear  a  uniform  temperature  for  all  look  down  angles. 

To  find  those  look  down  angles  which  give  significant  atmospheric  dependent  irradiance  values,  we 
can  choose  a  value  of  2  km  as  the  critical  path  length  such  that  in  a  maritime  midiattitude  summer 
environment,  the  atmospheric  transmission  is  between  0.5  and  0.6.  A  plot  of  look  down  angle  as  a 
function  of  vertical  height  can  be  determined  and  is  shown  in  figure  12.  Values  above  the  curve 
indicate  vertical  heights  and  associated  look  down  angles  where  the  airlight  can  contribute 
significantly  to  the  irradiance  at  the  observer  for  both  the  8  to  14  pm  and  3  to  4  pm  atmospheric 
windows. 


Figure  10.  Geometry  of  observer  relative  to  the  Earth  horizon  coordinate  system 


The  transmission  figures  were  obtained  from  Lowtran  6  for  the  8  to  12  pm  waveband. 


e  11.  L2as  a  function  of  viewing  angle 
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VIEWING  ANGLE  (DEGREES) 


Figure  12.  Vertical  height  as  a  function  of  viewing  angle 


5.  CONCLUSION 

The  statistical  infrared  signature  model  for  both  a  calm  and  wind  roughened  sea  enables  a  relatively 
simple  and  fast  numerical  technique  to  calculate  sea  radiance  values  to  produce  data  which  can  be  further 
used  in  processing  routines  incorporated  within  weapon  or  surveillance  system  algorithms.  A  further 
comparison  of  calculated  and  experimental  measured  sea  irradiances  for  different  wavebands  and  sea 
states  needs  to  be  made  for  a  higher  degree  of  validity  between  the  model  and  a  true  maritime 
environment. 

Using  Quadrature  was  one  way  of  numerically  evaluating  the  otherwise  awkward  integrals.  When  there 
is  a  high  degree  of  variation  in  the  sky  radiance,  such  as  broken  cloud  cover,  then  successively  higher 
point  Quadrature  must  be  used  if  the  same  numerical  accuracy  is  required. 


This  then  may  not  become  the  most  time  efficient  integral  solving  technique  and  an  alternative  integral 
evaluation  method  should  be  sought. 
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APPENDIX  I 

REFLECTED  ANGLE  DETERMINATION 

The  geometry  of  reflection  of  the  'sky'  off  a  wave  facet  to  the  observer  is  shown  in  figure  1.  a  is  the  angle 
made  between  the  observer  and  the  y  axis,  in  the  y  z  plane.  The  facet  has  a  maximum  slope  P  at  an 
azimuth  direction  a  taken  from  the  y  axis.  By  the  law  of  reciprocity,  the  reflected  ray  has  polar 
coordinates  v  and  p  where  the  angle  between  the  observer  and  reflected  sight  line  is  give  by  to. 

To  find  an  expression  for  the  component  of  the  maximum  slope  in  terms  of  maximum  slope  polar 
coordinates,  then  consider  figure  1.1  which  represents  a  cartesian  surface.  The  equation  of  the  tangent 
plane  at  a  point  'a'  on  the  surface  having  coordinates  (Xq,  Y0, 2 Q)  is  given  by 


Z 


Zq  +  Cx-Xq)  +(Y- Y0)f^r 


(1.1) 


where  the  normal  vector  at  'a'  is 


JH  -Jf 
dx  'dY  ' 


If  point  ’a’  becomes  the  coordinate  system's  origin,  then  equation  (1.1)  reduces  to: 


7  7  vaf  vaf 

Z  =  Zo+X37+YW 


(1.2) 


Since  grad  Z  gives  the  (x,  y)  coordinates  of  the  maximum  gradient  vectors  on  the  tangent  surface,  then  the 
above  equation  becomes: 


Z  = 


But 


Z  =  RtanP 


ERL-0473-TR 


24 


and  since 


-  ’  (SNSf 


R  =  tan  p 


Also,  any  vector  which  has  a  maximum  slope,  has  coordinates  given  by  x  =  Rsin  a  and  y  *  Rcos  a 
therefore 


sin  a  tan  p  = 


cos  a  tan  P  =  Zy 


where 


tan2  P  =  Z*2  +  Zy2 


thus  given  Zx  and  Zy,  p  and  a  can  be  determined. 

The  reflection  angles  v  and  ji  are  given  by  the  following  equations,  the  derivation  of  which  can  be  found  in 
reference  1. 


cosco  =  cos  p  sin  cos  a  sin  P  cos  $ 


cosn  =  2cos  P  cosw- sin  $ 


cotv  =  cota  +  ycscacscpsecwoos^ 


(1.3) 
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The  sign  of  the  second  term  in  equation  (1.3)  is  however  opposite  to  that  obtained  by  Cox  and  Munk  in 
reference  3. 


Figure  1.1  Wave  facet  as  tangent  plane 
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APPENDIX  II 

SEA  RADIANCE  INTEGRAL  EVALUATION 

Given  the  azimuth  and  elevation  angles  (a,$)  of  the  observer,  the  sea  radiance  integral  is  given  by 


=  -L  f 


tan  <(> 

|  (B  -  R[B  -  NS]}  cos  co  sec  | 


r  -(Zx2  +  zv2n 

P exp  —  JdZxdZy 


where  S  is  the  shadow  factor.  The  radiance  expression  in  the  above  integral  can  be  represented  by  a 
function  F(Zx,  Zy)  such  that: 


F(ZX,  Zy)  =  (B  -  R[B  -  NSI}  cos  co  sec  P 


so  the  integral  reduces  to: 


N(a,4»)  =  ^  Hi  +  U 


where 


r  ,  r  (zx2+zv2n 

1 1=//  FfZx,  Zy)  exp  -  -y-  dZxdZy 

-oo  -oo  L  Cr 


and 


r  ta,n  <t>  r  (zx2  +  zv2)i 

I2=  J  I  F(Zx,  Zy)  exp^-  -x  y--jdZxd Zy 


By  a  change  of  variables  the  integrals  Ij  and  I2  can  be  evaluated  using  Gauss  Quadrature  as  follows. 
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H.l  Change  of  variables  for  I 


letting 


Zx 

*  =  T 


Zy 
V  "  8 


I,  =  8 2  j  J  G(Y,V)e-Vdy  e'Y  dy 

-oo  -O 


where 


G(y,V)  =  F(8y,  5(-v»  eV-(-V)2 


n.2  Change  of  variables  for  I2 


by  letting 


Zx 

*  =  T 


Zy 

7— t. .  f 
Z  8  e 


where 
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I2  -  52  J  e  /  F(yS,  8(e  +  pe))  exp  -  (e  +  pe)2  dp  e'Y2  dy 
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